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: Abstract 

We present a systematic study of twist-4 light-cone distribution amplitudes of vector mesons 
in QCD, which is based on conformal expansion. The structure of meson mass corrections 
is studied in detail. A complete set of distribution amplitudes is constructed, which satisfies 
all (exact) equations of motion and constraints from conformal expansion. Nonperturbative 
input parameters are estimated from QCD sum rules. Our study suggests that meson mass 
corrections may present a dominant source of higher twist effects in exclusive processes. 
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1 Introduction 



The notion of distribution amplitudes refers to momentum fraction distributions of partons 
in a meson, in a particular Fock state, with a fixed number of constituents. For the minimal 
number of constituents, the distribution amplitude is related to the meson's Bethe-Salpeter 
wave function (pBS by 



The standard approach to distribution amplitudes, which is due to Brodsky and Lepage 0, 
considers the hadron's parton decomposition in the infinite momentum frame. A conceptually 
different, but mathematically equivalent formalism is the light-cone quantization ||. Either 
way, power suppressed contributions to exclusive processes in QCD, which are commonly 
referred to as higher twist corrections, are thought to originate from three different sources: 

• contributions of "bad" components in the wave function and in particular of components 
with "wrong" spin projection; 

• contributions of transverse motion of quarks (antiquarks) in the leading twist compo- 
nents, given for instance by integrals as above with additional factors of k\\ 

• contributions of higher Fock states with additional gluons and/or quark-antiquark pairs. 

In this paper we continue the systematic study of higher twist light-cone distribution 
amplitudes started in Ref. ||. In particular, we extend the analysis of || to include twist- 
4 distribution amplitudes and, most significantly, meson mass corrections. A preliminary 
account of some of our results has been reported in ||. 

Following ||, we define light-cone distribution amplitudes as meson-to- vacuum transition 
matrix elements of nonlocal gauge invariant light-cone operators. This formalism is perhaps 
less intuitive than the infinite momentum frame formulation, but it is more convenient for the 
study of higher twist distributions as it is Lorentz and gauge invariant. It allows all equations 
of motion to be solved explicitly, relating different higher twist distributions to one another. 
We will find that, much like in the twist-3 case ||, all dynamical degrees of freedom are those 
describing contributions of higher Fock states, while all other higher twist effects are given in 
terms of the latter without any free parameters. 

A systematic study of meson mass corrections presents the principal new contribution of 
this work. By counting dimensions, for any exclusive observable involving a large momentum 
transfer Q, power suppressed higher twist corrections have the generic structure 



twist-2, twist-3 and twist-4 operators, which have dimension 0, 1 and 2, respectively. The 
terms ~ m? do not involve any new dynamical information about the meson structure as com- 
pared to the leading twist terms, and are usually referred to as "kinematic" power corrections. 





-Um 2 • «0( 2 )» + m • «0( 3 >» + . 



Here m is the meson mass, ((O^)) and ((O^)) and ((O^)) are reduced matrix elements of 
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The structure of such kinematic corrections is well known for deep-inelastic lepton-hadron 
scattering, in which case they can be absorbed into a redefinition of the scaling variable || . 
The crucial observation leading to this "Nachtmann scaling" is that hadron mass corrections 
("target mass corrections" in this context) arise exclusively from the definition of the relevant 
leading twist matrix elements and do not involve new (higher twist) operators. This sim- 
plification does not hold in exclusive processes because there are additional contributions of 
operators containing total derivatives. Specifically, to twist-4 accuracy, in addition to Nacht- 
mann's corrections, there are also contributions of operators of type 



where is a leading twist operator. We find that contributions of the first type can be 
taken into account consistently for all moments, while contributions of the second type are 
more complicated and can be unravelled only order by order in the conformal expansion. 

The outline of this paper is as follows: definitions of and notations for distribution ampli- 
tudes are presented in Sec. 2 together with general remarks about specific features of the oper- 
ator product expansion (OPE) for exclusive processes and about conformal expansion. Section 
3 gives a general discussion of meson mass corrections for a simple example. The subsequent 
Sees. 4 and 5 contain a detailed derivation of chiral-even and chiral-odd distribution ampli- 
tudes, respectively. We take into account contributions of the leading and next-to-leading 
conformal spin and derive a self-consistent approximation for the distribution amplitudes, 
which respects the exact QCD equations of motion. The chiral-even and chiral-odd asymp- 
totic distribution amplitudes involve three nonperturbative parameters, and four additional 
parameters are required for the description of the leading corrections. The corresponding esti- 
mates are worked out using the QCD sum rule approach ||. On the basis of these estimates, 
we suggest that higher twist effects in exclusive processes are in many cases dominated by me- 
son mass corrections alone. The final Sec. 6 contains a summary and conclusions. The paper 
also contains three appendices in which we derive equations of motion for nonlocal operators, 
and derive and estimate the independent nonperturbative parameters that enter the twist-4 
distributions discussed here. 

Throughout this paper we denote the meson momentum by and introduce light-like 
vectors p and z such that 



The meson polarization vector e\t is decomposed into projections onto the two light-like 
vectors and the orthogonal plane as 



/il/i2-../in 



and 



2 pz 





(1.3) 
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Some useful scalar products are 



z ■ P = z ■ p = J(x- P) 2 — x 2 m 2 , 

p . e (A)_!i 2 . e (A) 

2pz 

e (A) . z = g (A) . x _ ^ 

We also need the projector onto the directions orthogonal to p and z: 

9^v = 9»v- —(p^+PuZf,), (1.5) 

and will use the notations 

a. = a M z M , a* = a^jf /{pz), (1.6) 

for an arbitrary Lorentz vector a M . 

We use the standard Bjorken-Drell convention for the metric and the Dirac matrices; 
in particular 75 = 27°7 1 7 2 7 3 , and the Levi-Civita tensor e^Ao- is defined as the totally anti- 
symmetric tensor with e i23 = 1- The covariant derivative is defined as D p = D ^ = d^ — igA^, 
and we also use the notation D ^ = d M + igA p in later sections. The dual gluon field strength 
tensor is defined as G^ u = ^e^ upa G pa . 



2 General Framework 

Amplitudes of light-cone-dominated processes involving vector mesons can be expressed in 
terms of matrix elements of gauge invariant nonlocal operators sandwiched between the vac- 
uum and the vector meson state, e.g. a matrix element over a two-particle operator, 

(0\u(x)T[x, -x)d(-x)\p'(P, A)), (2.1) 

where T is a generic Dirac matrix structure and we use the notation [x, y] for the path-ordered 
gauge factor along the straight line connecting the points x and y: 



[x, y] = P exp 



ig dt(x-y)„A fM (tx+(l-t)y) 







(2.2) 



To simplify the notation, we will explicitly consider charged p mesons; the distribution ampli- 
tudes of p° can be obtained by choosing appropriate isospin currents. 

The asymptotic expansion of exclusive amplitudes in powers of large momentum transfer 
corresponds to the expansion of amplitudes like ( |2.1| ) in powers of the deviation from the light- 
cone x 2 = 0. As always in a quantum field theory, such an expansion generates divergences 
and has to be understood as an OPE in terms of renormalized light-cone nonlocal operators 
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whose matrix elements define meson distribution amplitudes of increasing twist. To leading 
logarithmic accuracy, the coefficient functions are just taken at tree-level, and the distributions 
have to be evaluated at the scale [x 2 ~ x~ 2 . In this section we present the necessary expansions 
and introduce a complete set of meson distribution amplitudes to twist-4 accuracy. This set 
is, in fact, overcomplete, and different distributions are related to one another via the QCD 
equations of motion, as detailed in later sections. 



2.1 Chiral-Even Distribution Amplitudes 

We start with the matrix elements involving an odd number of 7 matrices, which we refer to as 
chiral-even in what follows. For the vector and axial vector operators the light-cone expansion 
to twist-4 accuracy reads: 

f P ( A ) T r 1 r m 2 x 2 

(0|«(x) 7 „d(-x)|p-(P,A)> = fp m p[ 1 ^ P nJ o due* p *^n(u) + -^A(ii) 



+ eW - R 



e w x 

Tx~ 



^ jf due* Px B(u) 



(0\u(x)^ 5 d(-x)\p-(P,X)) 



\ (7, - fj^±^) m p e^P a x p fdue^u). 
A V m p J Jo ^2.4) 



Notice that in order to calculate exclusive amplitudes to 0(1/Q 2 ) accuracy, terms of 0(x 2 ) 
have to be kept in the vector matrix element, but can be neglected in the axial vector one. For 
brevity, here and below we do not show gauge factors between the quark and the antiquark 
fields; we also use the short-hand notation 

£ = 2«-l. 

The vector and tensor decay constants f p and fj are defined, as usual, as 

(0|«(0) 7M d(0)|p-(P,A)> = f p m p e£\ (2.5) 

(0|u(0)(vd(0)|Ar(P,A)) = </J(eWp„-4 A) ^)- ( 2 - 6 ) 
The coupling fj is scale-dependent, with 

^ 2 »^> 2 »(H§) C> "' (27) 

with the standard notation C F = (N 2 - 1)/(2JV C ) and b = (11N C - 2n f )/3. 
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The expansions in ( |2.3| ), ( |2.4| ) involve several Lorentz invariant amplitudes, which we have 
to interpret in terms of meson parton distributions. Definitions of the latter involve nonlocal 
operators at strictly light-like separations and can most conveniently be written using the 
light-cone variables ( |1.2| ), and for longitudinal and transverse meson polarizations separately. 

Following |J, we define chiral-even two-particle distribution amplitudes of the p meson as 



\u{z)^^d{-z)\p (P,A)) = f p m f 



+ e 



(A) 



PZ Jo 

e^z 



due^ pz (P\\(u,p 2 ) 
m 2 [ due i ^ pz g 3 ( y u,p 2 ) 



(2.8) 



\u{z) llllh d{-z)\p-{P,X)) 



1 ( , _ f T m u + m d 

2V P Jp m n 



(2.9) 



The distribution amplitude <j>\\ is of twist-2, and g^' of twist-3 and g 3 of twist-4. All four 
functions <fi = {<f)\\, g±, g±\ g 3 } are normalized as 



Jo 



du(b(u) = 1, 



(2.10) 



which can be checked by comparing the two sides of the defining equations in the limit z p — > 
and using the equations of motion. 

Comparing fl2.8|) , Q2.9|) with the light-cone expansions in (|2.3| ), (|2.4j ), we easily find 

M(u) =g { l\u), 

C(u)=g 3 (u) + <P\\{u)-2g { l\u) 1 



B(u)=g^(u) 



(2.11) 



which is nothing but the tree-level OPE of the invariant amplitudes B, C, D in terms of meson 
distribution amplitudes. We will find (see also [[|) that all higher twist two-particle distribu- 
tion amplitudes do not present genuine independent degrees of freedom, but can be expressed 
in terms of three-particle distribution amplitudes. The same analysis will allow us to cal- 
culate the remaining invariant amplitude A, which accounts for the transverse momentum 
distribution in the valence component of the wave function. 

Three-particle chiral-even distributions are rather numerous and can be defined by the 
following matrix elements: 



{Q\u{z)gG lxv {vz)j a ^d{-z)\p (P,X)) = f P m p p a \p v e { ^ ~ Pn e ±l]A(v,pz) 
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■ _L_L $ $ 

■ * ■ $ ^ 


J_±.± T x (4) 
■ J_±± T 2 (4) 
• * • ± Tt ] 


■ ± S S 



Table 1: Identification of three-particle distribution amplitudes with projections onto different 
light-cone components of the nonlocal operators. For example, • _L_L refers to V'GaTxTsV^ 



+ f P m l -\Pii9tv ~ Pv9ap}®( V iP Z ) 



pz 



+ fp m lj^P4P^ -PvZ^(v,pz), (2.12) 



(0\u{z)gG flu (vz)i^ a d{-z)\p (P)) = fpm p p a \p v e { ^-p^eY v ]V{v,pz) 

,(A). 



+ fp m l [Pp9au ~ Pu9 a ^(v,pz) 

pz 



'(pz) 

where 



+ //^-^^PatP^-P^]*^^), ( 2 - 13 ) 



= y Pae-^ (a "- ad+ ^ 9 U(a), (2.14) 

etc., and a is the set of three momentum fractions a = {ad, ot u , a g }. The integration measure 
is defined as 

y Va = J dad J da u J da g 5 (l — . (2-15) 

The distribution amplitudes V and .4 are of twist-3, while the rest is of twist-4; we have not 
shown further Lorentz structures corresponding to twist-5 contributions^. Different distri- 
bution amplitudes can be separated by projecting onto particular light-cone components, as 
summarized in Table 1. 

For completeness, let us mention that also four-particle twist-4 distribution amplitudes 
exist, corresponding to contributions of Fock states with two gluons or an additional qq pair, 



■"■Note that we use a different normalization of three-particle twist-3 distributions compared to 
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of type 



^7.(75)^7.(75)^, ip G._l G._l7.^- 



Such distributions will not be considered in this paper for two reasons: first, it is well known || 
that four-particle twist-4 operators do not allow the factorization of vacuum condensates such 
as (ipip), (G 2 ). Because of this, their matrix elements cannot be estimated reliably by existing 
methods (e.g. QCD sum rules), although they are generally expected to be small. Second, and 
more importantly, the four-particle distributions decouple from the QCD equations of motion 
in the two lowest conformal partial waves. To this accuracy, therefore, it is consistent to put 
them to zero. Vice versa, nonvanishing four-particle distributions necessitate the inclusion of 
higher conformal spin corrections to distributions with less particles, which are beyond the 
approximation adopted in this paper. 

2.2 Chiral-Odd Distribution Amplitudes 

For chiral-odd operators involving a^ v and 1, the light-cone expansion to twist-4 accuracy 
reads: 



(0\u(x)a^d(-x)\p-(P,\)) 



WP M ) fdue* Px \<p 1 _{u) + 






(2.16) 




due l ^ Px B T (u). (2.17) 



The couplings f p and fj are defined in Q2.5| ) and (|2.6|). 

The corresponding distribution amplitudes on the light-cone are defined as 



(0\u(z)<r^d(-z)\p-(P,X)) 






i(fJ-f P 



m u + m d 




1 



(0\u(z)d(-z)\p~(P,\)) 



due^hfiu.p 2 ). 



m 



(2.19) 
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(s t) 

The distribution amplitude <J)± is of twist-2, h), ' are of twist-3 and h 3 is of twist-4. All four 
functions <p — {4 > ±^ ^4 j h 3 } are normalized as 



du(J)(u) = 1. 

Comparing ( |2.18| ) and (|2.19|) with the light-cone expansion ( ]2.16| ) and ( |2.17|) , we easily find 

,(*)/ 



1 1 

»t0) = h\* } (u) - -<p±{u) - -h 3 (u), 



C T (u) = h 3 (u) - (f)±{u). 



(2.20) 



As for chiral-even distribution amplitudes, only the twist-2 distribution 0j_ represents gen- 
uinely independent degrees of freedom, the others can be expressed in terms of three-particle 
distribution amplitudes. 

The three-particle distribution amplitudes are even more numerous than in the chiral-even 
case and can be defined as: 

{Q\u{z)a aP gG txl/ {vz)d{-z)\p'(P, A)) = 

■ rT 2 r (A) ± (A) J_ (A) J_ , (A) J_ im(4)/ \ 

, f T 2r (A) J_ (A) 1 (A) J_ (A) ±1^,(4), \ 

+ f P m p[Ppe^g l3u - Ppe^g au - p^g^ + p^g^T^ >{v,pz) 

fo m l 



+ 



' p p i 

pz 



\PaPpe±lzu - PpPp^ z v ~ PaPue^Zf, + p p p v e^ 3 a z ^T^ 1 (v,pz) 



+ 1 ^- R \PaPpe ( llzp - ppp^lza - p a p v e^z p + p0 u e^z a }T^'(v,pz), (2.21 



,(4), 



pz 



\u(z)gG fMU (vz)d(-z)\p (P,A)> = if^m 2 p [e { ^ v - e^p^S^.pz) 



\u{z)igG^(vz)-f 5 d{-z)\p (P,A)) = if^rri^^pu - e^lpy\S{v,pz). 



(2.22) 



Of these seven amplitudes, T is of twist-3 and the other six of twist-4; higher twist terms are 
suppressed. The relation of these distribution amplitudes to specific light-cone projections of 
the matrix elements was made explicit in Table 1. 

Also in this case there exist four-particle twist-4 distribution amplitudes which we do not 
consider for the reasons mentioned at the end of Sec. 2.1. 
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2.3 Conformal Partial Wave Expansion 



Conformal partial wave expansion in QCD [[], |It], [II], [12], |T^| parallels the partial wave ex- 
pansion of wave functions in standard quantum mechanics, which allows the separation of 
the dependence on angular coordinates from that on radial ones. The basic idea is to write 
distribution amplitudes as a sum of contributions from different conformal spins. For a given 
spin, the dependence on the momentum fractions is fixed by the symmetry. To specify the 
function, one has to fix the coefficients in this expansion at some scale; conformal invariance 
of the QCD Lagrangian then guarantees that there is no mixing between contributions of dif- 
ferent spin to leading logarithmic accuracy. For leading twist distributions the mixing matrix 
becomes diagonal in the conformal basis and the anomalous dimensions are ordered with spin. 
Thus, only the first few "harmonics" contribute at sufficiently large scales (for sufficiently hard 
processes). 

For higher twist distributions, the use of the conformal basis offers the crucial advantage 
of "diagonalizing" the equations of motion: since conformal transformations commute with 
the QCD equations of motion, the corresponding constraints can be solved order by order in 
the conformal expansion. Note that relations between different distributions obtained in this 
way are exact: despite the fact that conformal symmetry is broken by quantum corrections, 
equations of motion are not renormalized and remain the same as in free theory. 

The general procedure to construct the conformal expansion for arbitrary multi-particle 
distributions was developed in |QL O. To this end, each constituent field has to be decom- 



posed (using projection operators, if necessary) into components with fixed (Lorentz) spin 
projection onto the light-cone. 

Each such component corresponds to a so-called quasi-primary field in the language of 
conformal field theories, and has conformal spin 

j = ±(l + 8 ), (2.23) 

where / is the canonical dimension and s the (Lorentz) spin projection. In particular, / = 3/2 
for quarks and / = 2 for gluons. The quark field is decomposed as ip + = (l/2)^ip and 
■0_ = (l/2)^ip with spin projections s = +1/2 and s = —1/2, respectively. For the gluon 
field strength there are three possibilities: G,± corresponds to s = +1, G*± to s = —1 and 
both G±± and G* correspond to s — 0. 

Multi-particle states built of quasi-primary fields can be expanded in irreducible unitary 
representations with increasing conformal spin. An explicit expression for the distribution 
amplitude of a multi-particle state with the lowest conformal spin j = ji + ■ ■ ■ + j m built of m 
primary fields with spins jk is 

/ / \ T 2ji + . . . + 2j m 271-1 2jr' 2 -l 2j m -l /q Q/1\ 

<f) as {a 1 ,a 2 ,...,a m ) = a/ 1 a/ 2 . . . a^ m . (2.24) 

Here are the corresponding momentum fractions. This state is nondegenerate and cannot 
mix with other states because of conformal symmetry. Multi-particle irreducible representa- 
tions with higher spin j + n, n — 1, 2, . . ., are given by polynomials of m variables (with the 
constraint Yuk=i ak = 1 )' which are orthogonal over the weight-function ( |2.24| ). 
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3 Meson Mass Corrections 



The structure of meson mass corrections in exclusive processes is in general more complicated 
than that of target mass corrections in deep inelastic scattering where they can be resummed 
using the Nachtmann variable ||. For illustration, consider the simplest matrix element 

(0\u{x)^d{-x)\p~(P,X)) = 



f p m p (e w x) I due i{2u - 1)Px ^(u) + + 0(x 4 ) 



(3.1) 



We assume that x 2 <C A Q ^ D , but nonzero, (f>\\(u) is the twist-2 chiral-even distribution am- 
plitude and $(«) = mi[A(u) + (1/2) f dv f^dw C(w)} (c.f. ( |2.3| )) describes higher twist cor- 
rections, the "kinematic" contributions to which, due to the massive p meson, we want to 
calculate. 

Experience with inclusive distributions tells us that meson mass corrections are related to 
contributions of leading twist operators. Indeed, the conditions of symmetry and zero traces 
for twist-2 local operators imply 



(0| ujt{i Dx) n d \p~{P,X)) 



tw.2 



f p mM X) x) 



\pxY -^^il(P,)- 



((On)), 



(3.2) 



where [. . .] tw .2 denotes taking the leading twist part (subtraction of traces, in this case) and 
((O n )) is the reduced matrix element related to the n-th moment of the leading twist distribu- 
tion 



M^ = / du(2u-l) n <P\\(u) = ((O n )). 



(3.3) 



we find that the same 



Expanding fl3.1p at short distances x — > and comparing it with 
reduced matrix element gives a contribution to the twist-4 distribution amplitude: 



AfW = / du(2u- l) n $(u) 
Jo 



1 



n + 3 



m 2 p ((O n+2 )), 



(3.4) 



which is the direct analogue of Nachtmann's correction in deep inelastic scattering. 

As pointed out in [14], there exists an alternative possibility to describe the mass corrections 
by modification of the exponential factor in ( |3.1| ) rather than a contribution to the twist-4 
distribution amplitude. To this end, we write 



(0| u(x)td(-x) \p-(P,X)) = f p m p 



J tw.2 



du 



tw.2 



(3.5) 



where [. . .] tw .2 on the left-hand side correspond by definition to a subtraction of traces in all 
local operators in the Taylor expansion of the nonlocal operator at short distances. As shown 
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I5|, this definition implies that the nonlocal operator satisfies the homogeneous Laplace 

d 2 



equation 



dxridx* 1 



u(x)^d(—x) 



0. 



tw.2 



and the same condition has to be fulfilled by the function 



x)e 



iiPx 



(3.6) 
in order that 



tw.2 



Eq. ( p.5|) be satisfied. The solution can easily be constructed order by order in the {m 2 x 2 ) k 



expansion 14|. To twist-4 accuracy, we obtain 



V>x)e* Px 



(e w x) 



J tw.2 



e^ Px + 



m 2 x 2 ^ 2 



dtt 2 e^ Px + 0(x i ) 



(3.7) 



By taking moments, it is easy to check that Eq. (|3.1|) , with the higher twist distribution 
function <£>(«) given in (|3.4j), is equivalent to Eq. ( |3.5| ) with the substitution (|3.7|). 

The result in ( |3.4| ) is, however, incomplete. The reason is that in exclusive processes one 
has to take into account higher twist operators containing total derivatives, and vacuum-to- 
meson matrix elements of such operators reduce, in certain cases, to powers of the meson mass 



times reduced matrix elements of leading twist operators. In the present case, write [[15 

„2 rl 



u(x)qtd(—x) = u(x)$d(—x 



x 
+ — 

J tw.2 4 



dt 



2 



o 



dx a dx a 



u(tx)$d(-tx) + 0(x 4 



u(x)$d(—x) 



J tw.2 



X 

T 



dtr d 2 [u(tx)tfd(-tx)] 



+ contributions of operators with gluons + 0(x 4 



(3.8) 



where we used Eq. ( |A.9| ) to obtain the last line. In the matrix element we can make the 
substitution d 2 — > — m 2 . Expanding, again, at short distances, and comparing with the short- 
distance expansion of (|3.1|), we get an additional contribution to M„ , so that the corrected 
version of ( |3.4j) becomes 



Mi*) 



n + 3 



m 2 [{(O n+2 )) + «O n »] + gluons. 



(3.9) 



Assuming the asymptotic form of the leading twist distribution amplitude <p, <j>{u) = Qu(l — u), 
so that ((O n )) = 3/[2(n + l)(n + 3)], this equation for moments is easily solved and gives 



$(w) 



30m 2 (1 -u) 2 



^2 ^ 2 j- 

~m p + -m p C 4 



(3.10) 



where we have included the "genuine" twist-4 correction (term in £4) due to the twist-4 quark- 
gluon operator, see definition in Eq. (|4.6j ). The QCD sum rule estimate is £4 ~ 0.15 [fL6[], so 
that the meson mass effect on the twist-4 distribution function is by a factor 2 larger than 
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the "genuine" twist-4 correction. This is an important difference to deep inelastic scattering, 
where the target mass corrections are small. 

The above discussion is still oversimplified and does not provide us with a complete sep- 
aration of meson mass effects. The major complication arises because of contributions of 
operators of the type 

wy lix (iD IMl )...(iD» n )d . (3.11) 



tw.2 



Such operators can be expressed in terms of operators with extra gluon fields, which means, 
conversely, that certain combinations of qGq operators reduce to divergences of leading twist 
operators and give rise to extra meson mass correction terms. The corresponding corrections to 
twist-4 distributions involve, however, higher order contributions in the conformal expansion 
of the distribution amplitudes of leading twist and do not affect the result in ( ft.lOQ , which 
only includes leading conformal spinF]. A calculation of the next-to-leading corrections will be 
presented below. 



4 Chiral-Even Distribution Amplitudes 

In this section we derive explicit expressions for chiral-even distribution amplitudes of twist-4 
including the leading and next-to-leading contributions in the conformal expansion. We first 
give a short summary of the relevant results of || to twist-3 accuracy. We then discuss the 
conformal expansion of twist-4 three-particle distribution amplitudes and relate the coefficients 
to matrix elements of local operators. Only a few operators turn out to be independent, so 
that the number of nonperturbative parameters is reduced considerably. Finally, we calculate 
the twist-4 two-particle distribution amplitudes from the equations of motion (EOM). The 
quark mass corrections will be neglected throughout this section. 

4.1 Twist-3 Distributions 

A comprehensive study of p meson distribution amplitudes to twist-3 accuracy was carried out 
in |T7], |TB|, P|, and we begin this section by quoting the results relevant to the present paper. 

The leading twist- 2 distribution amplitude for the longitudinally polarized p mesons, <f>», 
is expanded as |T7|, |TB| 

<f>\\(u) = Quu 

The parameter a| is defined by the matrix element of a twist-2 conformal operator with 
conformal spin 3: 

(0\u^zdz) 2 d-hzdz) 2 uMP'(PA)) = ^(e^z)(pz) 2 f p m p al (4.2) 
5 35 

2 The reason why leading conformal spin is not affected is that the divergence of a conformal operator 
vanishes in free theory. 



i + 4 2^ 2 



(4.1) 
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and is scale-dependent: 



kg 2 ) = L^aiGu 2 ), 4 



25 



c, 



(4.3) 



where L = a s (Q 2 )/a s (p 2 ) and C F = (N 2 - 1)/(2JV C ), b = (11N C - 2n f )/3. The parameter a\ 
has been estimated from QCD sum rules; its value at the reference scale p = 1 GeV is given 
in Table g. 

The three-particle distributions of twist-3 read |L7], 



V(a) = 540 (3ul(a d - a u )a d a u al 



A(a) = 360 ( 3 a d a u a 2 



l+utf-(7a,-3) 



(4.4) 
(4.5) 



The dimensionless coupling ( 3 is defined by the matrix element 



(0\ugG^j a j 5 d\p (P,X)) = f p m p ( 3 



1 

3* 



+ ^f p m 3 p C4 



(4.6) 



where £4 is a matrix element of twist-4, which we will need below, while and uj a are defined 

as 



\ui{gG av z a [i Dz)-{iD z)gG au z a )d\p-(P, A)) = i{pzfe ( ^m p f p - ( 3 u\ + 0{z v ) (4.7) 

Zo 



and 



(0|^ 7 5 [iDz,gG^\ d- j(idz)u^ 5 gG^d\p-(P,X)) = -{pzf^ v m p f p ^C^f + 0{ 

respectively, where [ , ] stands for the commutator. 

The scale-dependence of the twist-3 parameters is given by (with Ca = N c ): 



UQ 2 ) = L^ /h Up 2 ), 7 3 C 



1 



Cf + 3Ca, 



(4.9) 



and 



U% - UJ 3 



a\Q 2 



V A\ M 
^3 - Ujf X 



pa 
J- q 



~ % A 3 F _ 3 A 
5 /m 4 /-^ 1 /"( I /nf 



(4.10) 



3 Note that we use a normalization of distribution amplitudes different from that in [jl7|, ||] . In the notation 
of Ref. | , w£ = w£ , C 3 A = Ca , and C 3 y = (3/28)0,^ . 
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Numerical estimates are given in Table |2|. 

Finally, the two-particle distributions of twist-3 are determined from the EOM 



g±\u) = 6uu 



^ ) H = ^(i+a + (^;+5c 3 ) (3^-1) 



+ 



9 

112 



°* 64 H 3 ~ U: 



(3 - 30£ 2 + 35£ 4 ) 



(4.11) 



4.2 Twist-4 Distributions 



Due to the odd G-parity of the operator in Q2.13|) , the distribution amplitudes $ and \& 
are antisymmetric under the exchange of and a u , whereas $ and ^ are symmetric. The 
distributions \& correspond to the light-cone projection J.G.* (see Table 1) and have the 
conformal expansion 

$(gt) = 120a u a d a ff ^ 00 + V^io(3a 9 - 1) + . . 



= 120a u a d a 9 + ^i (ad - « u ) + 



(4.12) 



respectively. Note that the leading spin contribution to \I> vanishes because of G-parity (for 
massless quarks). 

In turn, the distribution amplitudes $, $ correspond to the 7^G.j_ light-cone component, 
and before a conformal expansion can be performed, we first have to separate the different 
quark spin projections. To this end, we define auxiliary amplitudes: 



(Q\u{z)gG^(vz)wcT(ir{*d(-z)\p) = f P m 3 p — {p^L ~ Pv9^) $ n (y,pz) 



pz 

{0\u(z)gG^(vz)^ a ^.d(-z)\p) = f p m 3 p — (p^ v - p„g afl ) ®^(v,pz). 

pz 

The relation of $^ to the original amplitudes is given by: 

$(a) = - [$U + $JLT] («), 



(4.13) 



$(a) = - [$ n - $ iT ] (a), 
and their conformal expansion goes in terms of Appell polynomials: 



a, 



0oo + 0oi (a g - 3a rf ) + 0io ( a g - -a u 



(4.14) 
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$U(a) = QOa d a„ 



>')o + 0oi (a g - 3a u ) + 0i O ( a g - -a d 



(4.15) 



9 

T I IT 

where we have taken into account the symmetry properties, i.e. O q = 0g O , e ^ c - Combining 
everything, we obtain 



(1 - ag) + 0oi a g (l - ag) - Qa u a d 



+ 



a g {l-a g ) - -(a u + a a 



$(a) = 30a 2 (a u - a d ) 



0oo + 4>oia g + -0io(5a 9 - 3) 



(4.16) 



At this point, the expansion involves two parameters, 0oo and tpoo, to leading conformal 
twist accuracy, and four more (0io, 0oi, "0io, "0io) for the corrections. Our next task will be 
to relate them to matrix elements of local operators and find out how many coefficients are 
independent. 

For leading spin the answer is easily obtained by taking the relevant light-cone projections 
of the matrix element in 



nil 



2 1 

3C3+3C4, 



^00 



(4.17) 



Note that the "twist-4" distribution amplitudes receive contributions of both twist-3 and twist- 
4 operators. This is due to the fact that the standard counting of twist in terms of "good" 



and "bad" components as introduced in |20| is at variance with the definition of twist as spin 
minus dimension of an operator. See also the discussion in Sec. 2.2. of Ref. [|J. The parameter 



(4 in (|4.6|) can be explicitly defined as the matrix element of a pure twist-4 operator: 



(0\ugG, ulul5 d\p-(P,X)) = f p m 3 p e^C 



(4.18) 



Its scale-dependence is given by 



C 4 (Q 2 ) = ^ C/ W), 7 4 c = !|cV, 



(4.19) 



and the numerical value was estimated in |I(J from QCD sum rules, see Table g and App. C. 

The calculation of the next-to-leading order spin corrections is more involved and is pre- 
sented in detail in App. B. The main observation is that the four coefficients 0io, 0oi, "0io, "0io 
in fact involve only one new nonperturbative parameter, in addition to the ones defined above. 
We find: 



^01 



— 4 - — c 3 + — C3(^ + vY) + -C A u£. 
12 2 12 16 3 3 2 
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fp [MeV] 


<4 C3 




C 4 




198 ±7 


0.18 ±0.10 0.032 ±0.010 


-2.1 ±1.0 3.8 ±1.8 


0.15 ±0.10 


0.8 ±0.8 



Table 2: Parameters of chiral-even distribution amplitudes. Renormalization scale is \i 
lGeV. 



— all - 
12 2 


h 4 


Cs± 


16 


1 I, 


7 


c 3 ± 


9 


-a, — 

4 2 


12 


16 



21 

V'w = -7^2 - Y2 ^ 3 + iq^ 3 u ^ ~ ' 

foo= %&-^Cavi + ^C*<4, ( 4 -20) 
3 16 4 

where the new parameter ujf is defined as 

_ r ~ 1 4 - ~ 

(0|tt 7 ? 7 5 rf - -(id^ugG^j^dlp (P, A)) + (/i«i/) = 

= 2f p m%ui ( c Wp„ + eWP„) • (4.21) 

u;;f is estimated from QCD sum rules in App. C, with the result given in Table ||]. The one-loop 
anomalous dimension of the operator on the left-hand side of Q4.21|) is not known. 

A few comments on the structure of ( 4.20 ) are in order. First, as already mentioned, twist- 



4 distribution amplitudes contain contributions of operators of twist-3. Note that the twist- 
4 chiral-even distributions considered here correspond to longitudinally polarized p mesons, 
while the twist-3 parts appearing in ( |4.20| ) formally correspond to transversely polarized 
mesons. The physical reason why distributions with different polarization appear to be related 
is Lorentz symmetry: a longitudinally polarized p meson can be made transversely polarized 
by going over to the meson rest frame, rotating the spin and boosting back to the infinite 
momentum frame. The spin rotation, however, is not a member of the collinear conformal 
group. Because of this, the conformal structure of twist-3 additions to twist-4 amplitudes is 
rather complicated and does not match the naive expansion, similar to the case considered in 
App. B of 0. Formally, this is yet another complication of having a nonzero meson mass. 

Secondly, we find a term in a| that corresponds to the next-to-leading correction in the con- 
formal expansion of the leading twist distribution amplitude. This contribution thus presents 
an additional meson mass correction and appears, in technical terms, through the operator 
identity (see App. B) relating the divergence of a two-particle conformal operator to operators 
involving gluon fields. In this respect distribution amplitudes in exclusive reactions are funda- 
mentally different from inclusive distributions, which involve only forward-scattering matrix 
elements, so that matrix elements of operators with total derivatives vanish. 
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Third, an inspection of the numerical size of the entries in Table |2| reveals that the coef- 
ficients in ( f4.20|) are grossly dominated by the term in u\, which is a genuine twist-4 effect. 
We do not see any physical reasons for this dominance, but, if correct, it suggests that the 
above-mentioned complications may have an only marginal effect on phenomenology. 

Finally, we have to specify the two-particle twist-4 distributions g% and A defined in Sec. 2. 
They are not independent, but can be expressed in terms of <3> and \1/ by using the EOM, see 
Eqs. ( |A.10|) , ( |A.12|) . To next-to-leading accuracy, we obtain: 



03 (u) 



2 || 40 20 
7 a 2 + yC3-yC4 



c 2 1/2 (0 



, 27 ii 

1 ~28 tt2 



+ ^ 3 -^c 3 K + 3^})q 1/2 (0, 



» > i 4 / 1 II 10 25 _ 
A(u) = 30itV { - ( 1 + — af + — C 3 + y C4 



(4.22) 



If 9 ii 1 , 3 

+ - — a" H C 3 + - G 

5 V 14 18 S 8 s 



3^3 -w 3 



c 2 5/2 (0 



+ 10 



-24-^3 + ^3^-42(4^ 



jT dv J\w^l + C l 2 /2 {i w ) - 3&,(1 - £*) In i 



1 + & 



(4.23) 



with £ w = 2w — 1. The double- integral can of course be taken analytically: 

1+& 



^" dv jy w {l + C 2 1/2 (U - 3^(1 - In ]- 



= —uu(2 + 13uu) + -m 3 (10 - 15w + 6u 2 ) lnu + -u 3 (10 - 15m + 6m 2 ) In it. (4.24) 
10 5 5 

The resulting functions g% (u) and A(u) are shown in Fig. 1 by solid lines. The dashed curves 
are obtained by omitting the next-to-leading spin corrections (which is the approximation 
adopted in ||, pljj), and the dash-dotted curves correspond to taking into account the meson 
mass corrections only and neglecting all twist-3 and twist-4 matrix elements. It is evident 
that the contributions from next-to-leading order conformal spin are small in both cases. The 
mass terms clearly dominate g% (u) and constitute about half of A(it). 

We stress that the given expressions are exact, provided the three-particle distributions 
are taken in the above approximation. This means, in particular, that flOp and JQ3|) 
reproduce the exact second moments of g% and d 2 /du 2 A, i.e. the normalization of A, but 
the fourth moment of #3 (second of A) also includes (uncalculated) contributions from even 
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0.2 0.4 0.6 0.8 1 
u 



0.2 0.4 0.6 0.8 1 
u 



Figure 1: Two-particle twist-4 chiral-even distribution amplitudes of the p meson: g 3 (a) and 
A (b). LO means neglecting contributions of higher conformal spin for twist-3 and twist-4 
operators and the mass terms correspond to pure meson mass corrections. 



higher conformal spin operators. We have checked that the second moments agree with those 
obtained from Taylor expanding ( |2.3|) . 

Note that g% corresponds to the spin projection s = —1/2 for both the quark and the 
antiquark, and thus has a conformal expansion in Gegenbauer polynomials C l ' 2 (2u — 1), cf. 

oo 
fe=2,4,- 

The coefficients g% and g% can be read off ( 4.22 ). The conformal expansion of A is not 
straightforward and contains for instance logarithms. 



5 Chiral-Odd Distribution Amplitudes 

The construction of twist-4 chiral-odd distribution amplitudes parallels that for chiral-even 
distributions in the previous section. We first recall the results for distribution amplitudes of 
twist-3. Next, we derive the conformal expansion of three-particle distribution amplitudes to 
next-to-leading order in conformal spin and relate the expansion coefficients to matrix elements 
of local operators. The two-particle distribution amplitudes are then obtained from the EOM 
as derived in App. A. 
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5.1 Twist-3 Distributions 

The leading twist- 2 distribution amplitude for transversely polarized p mesons, <f>±, is expanded 

3 



as 



17, 18 



4>_l(u) = 6uu 



(5.1) 



The parameter is defined by the matrix element of a twist-2 conformal operator with 
conformal spin 3, 



(0\r lf r ± .{i D z) 2 d-^dz) 2 ua L .d\p-(P,X)) = ^ e ^(pz)ifJ a^, 



12 

35 



and is scale-dependent: 



4(q 2 ) = l^VGu 2 ) 



72 



10 r 

If 



(5.2) 



(5.3) 



The numerical value of has been estimated from QCD sum rules and is given in Table [5] 
(at the reference scale p = 1 GeV). 

The only existing three-particle distribution amplitude of twist-3, T, is given by: 



T(a) = 5A0( 3 ujJ(a d - a u )a d a u a 



a- 



(5.4) 



with cuj defined as 



(0\ua^(gG^Zp(iDz) - (iDz)gG^Zp)d\p~(P,\)) = ipz) 2 (e^z)m 2 p fj- (5.5) 

The parameter £ 3 was already defined in (|4.6| ). The scale-dependence of is given by || 

,t 25 



ul{Q 2 ) = L* l^l(p\ 73 - 



6 



C F -2C 



A- 



(5.6) 



Finally, we also quote the two-particle distributions of twist-3 as obtained from the EOM 



i + (^4 + ^3 T )(^ 2 -i: 



hf(u) = 3e + 3 -aie (5e - 3 ) + ^ 3 T (3 - m 2 + 350. 

Numerical values of the input parameters are collected in Table |3|. 



(5.7) 
(5.8) 



5.2 Twist-4 Distributions 

Due to the odd G-parity of the operator in ( |2.21|) , the distribution amplitudes are an- 
tisymmetric under the exchange of a u and a d , whereas S and 5* are symmetric. In order to 
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resolve the conformal structure of if, it is advantageous to exploit the fact that 0^75 is not 
independent of cr^, and to define a "dual" matrix element 



(0\u(z)ia a [3'y 5 gG ^ u (vz)d(— z) \p) = r.h.s. of (2.21) with T — > T. 



(5.9) 



One easily finds 



g~(3) = »y^3) y| 



(4) 



-T, 



3 > 



(4) 



n(4) 
■ 4 > 



?(4) 



(4) 



T, 



(-1) 



-T, 

i(4) 



('1) 



(5.10) 



We next note that the distributions and = — T 3 correspond to the Lorentz spin 
projection s = +1/2 for both quark fields and the spin projection zero for the gluon. Hence 

T[ A \q) = 120ti (a u - a d )a u a d a g , 

-120t w (a u - a d )a u a d a g . (5-H) 



T 3 (4, (2) 



For the distribution amplitudes S, S, if and T 4 , on the other hand, one has to separate 
different quark spin projections. To this end, we define auxiliary amplitudes 

(0\u(z) r ^gG^(vz)d(-z)\p-(P,X)) = i^m^e^ - e^p^(v,pz), 

(0\u(z) ri *i l5 gG, u (vz)d(-z)\p-(P,\)) = if T p m 2 p [e^ Pv -e^p,}S n (v,pz), (5.12) 
and, similarly, two more distributions and by replacing 7.7* — > 7*7.. The relations to 



the distribution amplitudes in (|2.21|), (|2.22|) are given by 



-if (a) 



S{a) 




+ S l \a)) 


S{a) 


= l(S U (a) 


+ S l \a)) 


T 4 (4) («) 


= l(S U (a) 


-S u (a)) 


T 4 (4) («) 


= i(^(a) 


-S iT (a)) 



(5.13) 



The auxiliary amplitudes are expanded in Appell polynomials as 



S (a) = 60a u a? 



■5oo + sio ( a g - - a u J + s 01 (a g - 3a d ) 



3a, 



soo + sio \ a g - - a d J + s 01 {a g 



(5.14) 
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and similarly for S*^ and . Here we made use of the symmetry of S under the exchange of 
the u and d quarks, i.e. s 00 ~ = SqJ, etc. 

From (|5.13|) and ( |5.14 ) it now follows immediately that 



o yen ) 


= 30a 2 . 


soo (1 


- a g ) 


+ SiQ 


S{a) 


= 30a 2 


SQO (1 


-a g ) 


+ s 10 


T 2 (4) («) 


= -30a 2 g (a u - 


ad) 


SQO + 




= 30a 2 g (a u - 


ad) 


soo + 



a r 



Qa u a d } 



ja 5 (l - a g ) - ^ (a 2 + a 2 ,) j + soi {a 9 (l - a 9 ) - 6a u a d } 



1 

2 

1 
2 



(5.15) 



At this point, the expansion involves two parameters of leading conformal spin, sqo and Soo, 
and six more (sio, Soi, Sio, Soij ^io, ^io) for the corrections. Our next task will be to relate them 
to matrix elements of local operators and find out how many coefficients are independent. 

For the leading spin, the answer is easily obtained by taking the local limit z —>■ of ( |2.22| ), 
so that 

(5.16) 



soo — Cli soo — (I 



with 



(0\ugG, u d\p-(P,X)) = if T p m%l{^ ] Pu 

(0|%G Mi/ i 75 d|p-(P,A)) = ifJm 2 p Cl(e^P v 
The parameters (J, (J renormalize multiplicatively P"Df : 



e (A) P ) 
ei X) P»). 



(5.17) 



(Cf + C 4 T ) (Q 2 ) = L 1+Ih (Cl + CI) (A 7+ = 3C A - I C F; 

(Cf - CT) (Q 2 ) = ^ 7_/6 (CJ - (?) (A 7- = 4C A - AC F . 



(5.18) 



The numerical values can be estimated from QCD sum rules, see Table f| and App. C. 

The calculation of the next-to-leading order spin corrections is involved and presented in 
detail in App. B. The main observation is that the six coefficients sio, soi, Sio, Soi, t w and tio 
involve three new nonperturbative parameters. We find: 



SlO 



SlO 



3 x 1 

— a 9 

22 2 8 



3 a 1 



c 3 ^+!<<Q (i) > 

28 



22^-8^ T -55«° (1, » 



7 

IT 

7 

IT 



<(Q (3) » 
«Q {3) » 



14 
~3~ 

14 



<(Q (5) », 
«Q (5) », 
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*oi = ^ 4 + 1 C*% + ^ ((Q {1) )) - y 2 ^ (3) )) + 1 ^ 



'10 



~ 4 + ~ C*>Z - ^ + y 2 ^ (3) )) + 1 ^ 

A 4 _ A CaW T _ JA //q(i)\\ + 1H (/qo))), 

44 2 16 3 220 vx 77 44 ^ 77 



(5.19) 



The above relations involve the three parameters ((Q^)), ((Q^)) arid ((Q^)), which can be 
defined as matrix elements of the following operators: 



-iu V a (cr^gGnp - a np gG ip )d + 7uV a g[G - ij 5 G] iv d - — d a ug[G - ij 5 G]^ v d, 



(3) 



1 _ 1 _ 

ztt cr vp g[G + ij 5 G) a p d - - uV a g[G + i"f 5 G]cv d ~ 3 uV v9[ G + HsGka d 



+ ^ d a ug[G + ijsG]^ d + ^ d v ug[G + ij 5 G] ia d \ - {£ <-> 77}, 



(5.20) 



where we used convenient short-hand notations (see |19| ) for the covariant derivatives: V a 

G^y = G pu D a — D a G pv acting on quark fields only, and V a G pu = [D a , G pu ] acting on gluon 
fields only. The reduced matrix elements ({Q^ 1 ')} of these operators are defined as 



<o|QS>-(P,A)> 



ef ] (p a P v ~ \ m 2 p g a A - e^ A) ( P a P^ - \ m 2 p g a ^ 



+ (e™g cm -eWg ae )((B®)), 



(5.21) 



where ({Q^ 1 ')} is of twist-4 and ((R®)) of twist-5. The operators renormalize multiplicatively 
and their one- loop anomalous dimensions are known [I9|; to obtain the scale-dependence of 
the matrix elements ((Q^)), one has to subtract the anomalous dimension of fj, Eq. ( |2.7| ), so 
that 

<(Q W »(Q 2 ) = L'^ /b «Q ( W) 



11 r, 10 n 

7q(i) = -4C F + — C A , 7q( 3 ) = — C F , 7q( 5 ) 



-C F + 5C A . (5.22) 
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/J [MeV] 4 C3 

160 ± 10 0.20 ±0.10 0.032 ±0.010 7.0 ± 7.0 



Table 3: Parameters of twist-2 and twist-3 chiral-odd distribution amplitudes. Renormaliza- 
tion scale is yu = 1 GeV. 



U 


a 




«Q (3) » 


«Q (5) » 


0.10 ±0.05 


-0.10 ±0.05 


-0.15 ±0.15 









Table 4: Parameters of twist-4 chiral-odd distribution amplitudes. Renormalization scale as 
above. 



Numerical estimates for these matrix elements were obtained in [[0| and App. C and are 
collected in Table £|. 

Finally, we have to specify the two-particle twist-4 distributions h 3 and A^ defined in 
Sec. 2. They are not independent, but can be expressed in terms of S, and by using 
the EOM, see Eqs. ( |A.13|) , ( |A.14|) . To next-to- leading accuracy, we obtain: 

3 i - ~ ' >t I 7t\ 1 n l / 2 t 



h 3 (u) = l + \ -i + ~ 4 - met +Q)} C£"(0 

+ \~4-^Cs"t\C 1 A0, (5.23) 

At(«) = so^v ^ (l + ~ 4 + y (I - f Ci) + (| 4 + ^ Cs^ T ) cTtt) 

^4-^ 3 ^ + ^((Q (1) )) + ^((Q (3) )) 

x (uu(2 + 13ww) + 2u 3 (10 - 15m + 6m 2 ) lnw± 2u 3 (10 - Ihu + 6m 2 ) In it) .(5.24) 

In Fig. 2, we plot h 3 and A^ as functions of u, showing full results and the contributions from 
leading order conformal spin and mass correction terms separately. Like in the chiral-even 
case, the mass terms dominate h 3 (u) and constitute approximately one half of At(u). 

We stress that the given expressions are exact provided the three-particle distributions 
are taken in the above approximation. This means, in particular, that ( |5.23|) and ( |5.24|) 



reproduce the exact second moments of h 3 and d 2 /du 2 AT, i.e. the normalization of At, but 
the fourth moment of h 3 (second of At) also includes (uncalculated) contributions from even 
higher conformal spin operators. We have checked that the second moments agree with those 
obtained from Taylor expanding ( |2.16 ). 
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Figure 2: Two-particle twist-4 chiral-odd distribution amplitudes of the p meson: h 3 (a) and 
At (b). LO means neglecting contributions of higher conformal spin for twist-3 and twist-4 
operators and the mass terms correspond to retaining meson mass corrections only. 



Note that, like #3, h 3 corresponds to the spin projection s = — 1/2 for both the quark and 
the antiquark, and thus has a conformal expansion in Gegenbauer polynomials C 1//2 (2w — 1), 
cf- CT): 

00 

h 3 (u,p 2 ) = l+ hf\f)Cl /2 {2u-l). 

fc=2,4,... 

The coefficients h% and can be read off ( p.23| ). The conformal expansion of At is more 
complicated. 



6 Summary and Conclusions 

In the present paper we have studied the twist-4 two- and three-particle distribution ampli- 
tudes of vector mesons in QCD and expressed them in a model-independent way by a minimal 
number of nonperturbative parameters. The work reported here is an extension of our earlier 
paper on twist-3 distribution amplitudes || . The one ingredient in the approach is the use of 
the QCD equations of motion, which allow us to reveal interrelations between different distri- 
bution amplitudes of a given twist and to obtain exact integral representations for distribution 
amplitudes that are not dynamically independent. The other ingredient is the use of confor- 
mal expansion: analogously to partial wave decomposition in quantum mechanics, it allows 
one to separate transverse and longitudinal variables in the wave function. The dependence 
on transverse coordinates is represented as scale-dependence of the relevant operators and is 
governed by renormalizat ion-group equations; the dependence on the longitudinal momentum 
fraction is described in terms of irreducible representations of the corresponding symmetry 
group, the collinear conformal group SL(2,R). The conformal partial wave expansion is ex- 
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plicitly consistent with the equations of motion since the latter are not renormalized. The 
expansion thus makes maximum use of the symmetry of the theory in order to simplify the 
dynamics, which is related, in the perturbative domain, to renormalization properties of the 
relevant operators. 

The analysis of twist-4 distribution amplitudes is complicated by the fact that the twist- 
4 terms are of different origin: there are, first, "intrinsic" twist-4 corrections from matrix 
elements of twist-4 operators. There are, second, admixtures of matrix elements of twist- 
3 operators, as the counting of twist in terms of "good" and "bad" projections on light- 
cone coordinates does not exactly match the definition of twist as "dimension minus spin" 
of an operator. There are, third, meson mass corrections, which one may term kinematical 
corrections, that come, on the one hand, from the subtraction of traces in the leading twist 
operators and, on the other hand, from higher twist operators containing total derivatives 
of twist-2 operators. Meson mass corrections of the first kind are formally analogous to 
Nachtmann corrections in inclusive processes, while the contribution of operators with total 
derivatives is a specific new feature of exclusive processes, which makes the structure of these 
corrections much more complex. 

Our final results are collected in Sees. 4 and 5. We present a complete set of distribution 
amplitudes that is consistent with the QCD equations of motion and has a minimum number 
of nonperturbative parameters whose numerical values are estimated from QCD sum rules. It 
turns out that the meson mass corrections are the dominant ones in all two-particle twist-4 
distributions, which is in contrast to what is observed in deep-inelastic scattering and welcome 
from the phenomenological point of view, as the higher twist matrix elements, £3,4 etc., come 
with considerable numerical uncertainties. 

The results of our study are immediately applicable — and, in fact, have already been 
applied [^TJ - to processes such as exclusive or radiative B decays and hard electroproduction 



of vector mesons at HERA. 
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Appendices 

A Equations of Motion 
A.l Operator Identities 

In this appendix we collect exact operator identities, which can be derived using the approach 
of [|lj| and which present a nonlocal equivalent to the equations of motion for Wilson local 
operators. The basic idea is to study the response of nonlocal operators to total translations 
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and/or the change of the interquark separation along the light-cone. For convenience, we work 
in the Fock-Schwinger gauge x^A^x) = 0, so that 

[x,-x] = l, A u (x) = / dvvx a G alx (vx). 

Jo 

All operator relations can be made manifestly gauge invariant by restoring the path-ordered 
gauge factor between field operators at different points in space-time. 
For the general two-particle operator, we can write 

d *+ f 1 

u{x)Td(—x) = — u(x)r Dfj, d(—x) — ij dvvu(x)x Q gG a ^(vx)Td(—x), (A.l) 



9x u j _i 



where Y is an arbitrary Dirac matrix and D^—Dfj, — D u — (d —iA(—x)) u — (d +iA(x)) u . The 
derivatives act on the arguments of the quark operators. 

In a similar way we can calculate the derivative with respect to the total translation: 

<9^{u(x)r<i(— x)} = u{x){Dn + Dp)Td(— x) — i J dvu(x)x a G au (vx)Td(—x), (A. 2) 



where, by definition, 

f) 

(A.3) 



d 

dn {u(x)Td(-x)} = —— {u(x + y)[x + y, -x + y]Td(-x + y)} 
9y u 



Here it is important to keep the gauge factors, which give a nonvanishing contribution: 

d u [x, —x] = iA u (x)[x, —x] — i[x, — x]A^(— x) — ij dv [x,vx]x a G afl (vx)[vx, —x]. 



For chiral-even operators, T = { / y fl , 7^75}, the first terms on the right-hand side of Eqs. ( A.l 
( |A.2| ) vanish by virtue of the massless Dirac equation, so that 

d f 1 

u(x)^ u (j 5 )d(-x) = - ij dvvu(x)x a gG aiJi (vx)^ u (^ 5 )d(-x), (A.4) 



dx u 



dp{u{x)^^)d(-x)} = - ij dvu(x)x a G a p l (vx)^ a ('j 5 )d(-x). (A.5) 
For chiral-odd operators, T = {1(75), 0^(75)}, on the other hand, we can use the identities 

a^Dp=i(Du + Du), (5„ + 5„)<V = i 5„, (A.6) 
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and get, combining Eqs. ( |A.1| ) and ( |A.2| ): 



d f 1 

d IJ ,u(x)a lxu (j 5 )d{-x) = -i — u(x)(j 5 )d(-x) + J dvvu(x)x p gG pv (vx)(j 5 )d(-x) 



i / dvu(x)x p gG m (vx)a flv ('j 5 )d(-x), (A.7) 



d _ _ f 1 

-Q^u(x)a^(^)d(-x) = -id u u(x)("f 5 )d(-x) + J dvu{x)x p gG pv {vx){^)d{-x) 



-i J dvvu(x)x p gG p ^(vx)a fJ , u ('j 5 )d(-x). (A.8) 

This method is general and can also be used for calculating the second derivative. In particular, 
the following formula is useful: 

d 2 



dXrydX C 



u(x)T d(—x) = —d u(x)T d(—x) + u(x) [TaG + crG T]d(—x) 



d f 1 f 1 

— 2ix u — — / dvv u{x)TG vp (vx)d(—x) — 2ix u d p I dv u{x)TG up {vx)d{—x) 
u x nJ-i J -l 

/l PV 
dv J dt(l + vt)u{x)Tx^x u G pp {vx)G p u {tx)d{-x) 



+ ix u J dv(l + v 2 ) u(x)T[D^ G>" v \{vx)d{-x), 
where [D^, G p u ] = —t A {ip^ v t A ip) , assuming summation over light flavours ip. 



(A.9) 



A. 2 Relations Between Distribution Amplitudes 

We are now in a position to derive relations between two- and three-particle amplitudes. 
Sandwiching ( |A.4|) between the vacuum and the p meson state, we find 



d f u f u 1 
g 3 (u) = 0|| (u) - 2 — / da d / da u — {2$(a) + *(a)} . 
du Jo Jo a g 

To arrive at Eq. ( |A.10| ), the following formula proves useful: 
J dvv k J Va J-'(a) exp[— ipx{a u — ad + va g }] = 



(A.10) 



ru ru 

du exp[i^px] / dad / da u 
Jo Jo 
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a„ 



(a d - a u - £) 



J 7 (a), (A.ll) 



which is valid for an arbitrary function J 7 (a). 

On the other hand, taking the matrix element of ( |A.5|) , and eliminating g 3 by virtue of 
( |A.10D , we obtain 



l£ s A(u) = A(gl(u)-<f>\\(u))+A-f f da d f da u — {2$(a) + 
o au au Jq Jq (Xg 

d 2 /•«/■«! 
+ — / da d da u — (a u u-a d u){2$(a) + ^(a)}. (A.12) 
du J J a 



can then be obtained as 



A(u) 



U PV 



JO 



d 2 



dv / dw -— A{w) 



dw 2 



Two more relations are derived in a similar manner between chiral-odd distribution am- 
plitudes using the operator identities in ( A.7 ) and ( |A.8|) . From (|A.7|) : 



h\*\u)--(h 3 (u) + ( f> ± (u)) 



du 



u ru 



JO 



da d I da u 



ad a 2 u t S (a)--(-lf\a I 
ai a, 1 



,(4), 



(A.13) 



and from (|A.8| ): 



U fV 



1 



4 / dv I dw[h 3 (w) - (j)±(w)\ - - A T (u) - / dv (2v - l)[<f>±(v) + h 3 (v)\ = 
Jo Jo 2 Jq 

= [ U da d [ U da u - ! ad - a "-Z ( T ?\a) - T 3 (4) (a)) - S(a)\ . (A.14) 
Jo Jo oc g y a g V / J 

B Short-Distance Expansion of Distribution Amplitu- 
des: Relation to Local Operators 

In this appendix we calculate the next-to-leading corrections to the conformal expansion of 
twist-4 three-particle distribution amplitudes as given in Eq. (|4.20| ) and ( |5.19|) . 



B.l Chiral-Even 

Our general strategy will be to consider matrix elements of the relevant operators with all 
Lorentz indices open. Taking different light-cone projections we will relate coefficients in the 
conformal expansion of distribution amplitudes to the invariant Lorentz structures and then 
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identify the relevant contractions of indices that relate 0io, <poi, V'io and ipio to matrix elements 
of independent twist-4 operators. 

To the next-to-leading conformal spin accuracy, we need local operators of dimension 6 with 
one quark- ant iquark pair, one gluon field and one additional covariant derivative. Taking 
into account G-parity, only two operators can contribute, apart from operators with total 
derivatives: 



(D 



u(i Dp gG^v + gGfjayi Dphalsd, 



°ap^ = u (~ Dp gG^ + gG^ Dp)j a d. 
For each of them we write down a general Lorentz decomposition: 

<0|O$>-(P,A)> = (eW 



PaPpPv - ^ m p (P a g Pv + Ppg au ) - - m p P v g a p 



(b.i; 



5 1 
P a PpP» - ^ m 2 p (P a g Pfl + Ppg ap ) - - m 2 p P p g a p 



1 2 

24 p 



e^\g P uP„ - g^Pu) + e^ A) (^P M - 9aii P v ) \ A®f, 



+ P a {efg Pv - ePgpJB® + Pp(e^g m - e^g^C® 
+ e^(P M ^ - P v g Pll )D® + 4 A) (P^- - P»g a ,)E {i) 



-g a p(P^-P v e^)F®. 



fB.2) 



Here is of twist-3 and can easily be related to an integral over the twist-3 distribution 
amplitudes A and V, respectively. Using (|4.7|) and (|4.8|), we find 



(pz) 3 e^m p f p AM = (0\O {1 {\p) 



7 2£ 



3 > 



(B.3) 



:b.4) 



To project onto the intrinsic twist-4 contributions, we must replace one "dot" projection 
by a "perp" projection in ( |B.3| ) and (ET1), which yields 



|o£» = f p m%ez)(pz)g 



±± \ g ( _ 3 ^ 3 + 3 _ 14 ^ 01 + ^ 10 ' 



(0|Og» = /X(ez)(pz)^ 



■j6 + gC) + n (*« 



(B.5) 
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whereas direct contraction of ( |B.2| ) gives 

(0\Of.. ± \ P ) = ~rn 3 p f p (ez)(pz)gi ± AW + {ez){ P z)g{ ± (C« + E®) . (B.6) 

Once CW 4- £;W are known in terms of £3, C4 an d ^4 > these two equations serve to determine 
0oi and 0io- 

In order to determine C® + we first have to introduce some more matrix elements: 
(0|Og ? » = (e p P v + e v Pp)X® + (e p P„ - e v P p )X® , 

(®\0% v \p) = (e a P v + e v P a )Y® + (e a P u - e v P a )Y®. (B.7) 
By construction, X + and Y + are of twist-4 and X_ and Y_ of twist-5; thus, to our accuracy: 

X® = Y® = 0. 
Also note that, by definition, Eq. ( f4.21|) : 

X? = U (W - ^) ■ (B.8) 

Now, by contracting ( |B.2| ) with g a p, etc., we find a set of linear equations relating A, B, ... to 
X + and Y + , which can be solved to give 

B = D= ±(X + -3Y + )f p m% 



C = E = -±(X + -3Y + )f p m% 



F = 0. (B.9) 

We recall that contributions from twist-5 operators are neglected in these solutions 

Of the remaining three unknowns 
by observing that D^G^ U = 0, so that 



Of the remaining three unknowns X^ , Y^ and Y^ , Y+ can be obtained rather easily 



(0\O { X + O^Jp) = P^(0HG^lal 5 + G ialul5 )d\p) 

= -{e u P a + e a P u )f p m 3 p Q ( 3 + ^ Ca) > 

which means 

Kf 1) = -^C3-^C4- (B.10) 
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(2) (2) 

In order to determine the remaining parameters X\_ and Yl , we make use of the operator 



identities 



and 



-^d p E pa p 



n (2) n (2) _ n (l) n (l) n (l) 



fB.ir 



-Uiujp |g p/3 D Q - D a Gpp + (a /3)| d - Ad p u(^ p G ap + ^ a G Pp )^d 



where 



8 ~ 8 ~ 28 ~ 

- g d(3U-f a G aa -f 5 d - - d a u-f a G a pj 5 d + — g a /3d p u-f a G ap d, 



E, 



fia/3 



15 «-» <-» 3 

Y«7m D a Dp d - -d a dpu-i p d - traces 



- symmetrized 

is a leading twist-2 conformal operator. Taking matrix elements, we find 



X 



(2) 



Y. 



(2) 



X 



(1) 



Y 



(i) 



+ ' 



(B.12) 



(B.13) 



and 



- m 2 p 4 = -24/ p mJxf + 4/„mJ ( r + - «, 



1 



3 /p™pC4 



(B.14) 



With X® and from Eqs. (gj), flBlOD and (gOg ), we get C + £ from flBj), and thus 
0io and 0oi from ( |B.5| ) and ( |B.6| ), see the first two lines of Eq. ( |4.20| ). 

Note that it is precisely operator relations of type (p,12j ), where the divergence of a leading 
twist conformal operator is expressed as a certain combination of quark-quark-gluon opera- 
tors, that make the analysis of meson mass corrections to twist-4 distribution amplitudes so 
complicated. This divergence vanishes in a free theory, as expected. 

The determination of the remaining parameters -010 an d V'io is now fairly easy: introducing 
a different "bad" component in (B.3|) and (B.4Q, we find 



(B.15) 



(0|OE|P) = f P m%ez)(pz) \ {a^ - (x« + if) 



On the other hand, taking proper integrals over distribution amplitudes: 



(0|O«|p> = (ez)(pz)f p ml ( ^oo - ^io 



(0|O.(. 2 )|p) 



21 



f p m 3 (ez)(pz)ip w . 



(B.16) 



By equating (B.15| ) and (B.16| ), we obtain the last two lines of Eq. (|4.20|) . 
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B.2 Chiral-Odd 



The calculation of next-to-leading order spin corrections to chiral-odd distribution amplitudes 
essentially parallels the calculation of similar corrections to the photon distribution amplitude 
in Ref. To follow this analogy, it is convenient to express the corrections in terms of 



3 1 

Vi = -^(s w + 2s 01 ), T] 2 = (sio -2s i) 



(B.17) 



and the corresponding "dual" quantities rji, rj 2) instead of the coefficients soi, sio>sio> «oi in 
the expansion over Appell polynomials. 

Expanding ( |2.21| ) and (|2.22| ) in powers of (pz) to first order, we obtain 



(0|tZ V. <T a pgG ia/ d\p~) = fjm 2 p — [p a p^g pv - ..]ipz[ -— Cs^J 
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12— 1 ( \ 1 3 

- [paP^Z v -•••]— ^lO + [VaP^ - . . .] — f - Cj + — Vl + ^4 7 

(0|M#£>.G^d|p~) = fjm 2 p pz(e^p u - e^p M ) {^4 - ^ (?7i + j 



0(twist-5), 
(B.18) 



(0|m^.G^rf|p-) = f p m^z{e^p v - e^ M ) - ^ (t/j [> • (B.19) 



Comparing these expressions with the reduced matrix elements of the conformal operators, 
Eq. ( |5.21| ), we immediately find 



((Q 



(5)\ 



j^(Vi + V2) - ^(Vi + V2), 



while contracting (|B.18| ) over g@ u , we find 



((Q (1) )) = ~ (ho - fio) - ^foi - vi) - \ (i» - fc). 



Next, we introduce one more operator 
Q (2) 

Using ( B.18|) and ( B.19|) , we easily find 



„., (/ - «M V<* {a^ p gG vp - a vp gG ip )d + ~ uV a (gG^ - igG^ v j 5 )d. 



«Q (2) » = ^ (*io - M - ^ (771 - >ji) - ^ (V2 - m). 



(B.20) 



(B.21) 



(B.22) 



^B.23) 
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is actually not independent, but related to via an important operator identity derived 
in pi: 4 ( ^ 

Qa\r, = 2 Qafo + 3Q ®p {Ol,mp ~ 0%£ ap ) , (B.24) 

where O 2 is the leading twist-2 conformal operator 

15 <-» <-> 3 
°a-.. = y V. V. d - - d 2 ua a .d. (B.25) 

The matrix element of the conformal operator on the right-hand side of Eq. ( p.24j) is equal to 



(0\O. 



2 ^ ap \p ) = if J (~y a 2^ ^PaP P (e^P v - e v P^) - i m 2 p g ap (e^P ri - e v P^ 



1 



+ ^ m P [P»(9^ P en ~ g VP et) + P P {g^e v - g m e^)} J , (B.26) 

so that 

((Q (2) )) = ^((Q (3) ))-^4- (B.27) 

At this point we have established three relations for the six independent parameters. Three 
more relations follow from the analysis of the most general matrix element 

{Q\Oa,va>£i\p~) = (0|w V Q (7tf»igGt v d\p~~) 

= P a [e^P^r, - P v g^) - e v (P^g pv - P T ,g P n)] A 
+ P a [P p (e^g uv - e v g u( .) - P v {e i g m - e v g^)] B 
+ e a [P^P^gvr, - P v 9u() ~ P v {Pt9 m ~ Pr,9i^)\ C 
+ {g ap Pv - g*vP P ){e£, P v ~ e v P d D + ig<Pf) - g an Pt){e p Pv - e u P^)E 
+ e a {g^gu V - g m gvi)F + [g a n(g&e v - fiv e ?) - g^ig^ - g m e^) G 

&v gr/u^fi) gar) {g^u - g&e-n)] H - (B.28) 

By projecting onto different light-cone variables, we find a set of linear equations for the 
coefficients A,. . . , E (F, G, H are of twist-5 and thus not relevant for the following discussion): 

A + B + C = ffml^&vl, 



4 Note that we obtain a different sign in front of the total derivative operator on the right-hand side as 
compared to Eq. (4.24) in Eg]. We thank G. Stoll for checking this equation. 
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B 



21 Jp p 



10: 



A = fJ m U^t + i;Vi + i : 7V2 



1 

42 



3 

14 



A-E 



1 l t 1 



5-^ = /p T ^(gC 4 T +^m + T7^)- 



3 

14 



(B.29) 



Contracting ( |B.28| ) with ^ a ^, we find 

{Q\uD v G^d- d v uG^d\p~) = (0\O PtPV ^ v \p~) 

= P u \Pte ri - P v es)(B - C - 3D - E) + [e^ - e v g^](m 2 p B + F-3G-H). (B.30) 
Using ( B.19Q , this can be translated into 



so that finally 



(*io - * 



io J 



C3W3 ~ 4 V2- 



The same analysis can be performed for the matrix element of the dual operator 

O a^jjbu^rj Va <7 P yG ^d , 



which results in the "dual" version of the relation ( |B.32| ): 

2 ~ 3 

o (*io - *io) = o Cs^J - 4r 72- 



fB.31} 



(B.32) 



^B.33) 



(B.34) 



To obtain the last relation from which 771,^2,^1,^2,^10,^10 can be extracted, we use the 
identity: 



(3) 



-- I (5) 
3 V ^ 2 



TjCt 



Ul 



V$ cr^g^G + i^G) ap - V v cr(. p g(G + ^G) 



OLp 



^B.35) 



where Ot = u[G + i^G\ n d. The matrix element of the first three terms on the right-hand 



side is: 



(OI-3C) 



~ (•••) + ~ (• ■ Ob") = ^« - esP v )fJm 2 p { «Q ( 



(C 4 T + C 4 T ) • (B.36) 
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For the remaining term on the right-hand side we use that 

1 



■ g a7) aG, 



and observe that the term in aG has zero matrix element over the p meson. Thus, using 
( jig ), we find 



(0\ui V s (cr vp gG ap + a ap gG vp )d -(£<-> rf)\p~) = P a (P^e v - P v e^)[2A + IB - AC - 3D - 3E). 

(B.37) 

This gives 



((Q 



(3)x 



2A + 2B - AC - 3D - 3E + ((Q (5) )) - i (Cj + CD, 



and with A, . . . , from (B.29| ), finally: 



(3) » + «Q (5) » = y a C^l + ~ (ho + fio) - ^ fa + 9^ 2 ) - i fa + 9?j 2 ) 



(B.38) 



(B.39) 



7 V " 1U ' " luy 14 v,x ' 14 

The eight relations (|B17l ), (|B~20|) , ( ggg) , ( |B~23D , (p~27| ), ( jBgg ), QBgg ), (|B~39|) yield s 10 , s i, 
sio, sbi, *io, *io as given in ( |5.19| ). 



C Numerical Estimates: QCD Sum Rules 

In this appendix we estimate the independent nonperturbative parameters from QCD sum 
rules. The sum rules for chiral-odd matrix elements can be adapted from the analysis of the 
photon distribution amplitudes in Ref. [19], while those for chiral-even matrix elements are 
partly available from Ref. |IB[ , partly new. The numerical results are collected in Tables [| [| 
and |. 



C.l Chiral-Even 



To leading conformal spin accuracy, we need the single parameter £4 ( |4.18| ). This matrix 
element was discussed at great length in fi~6fl . The best estimate comes from considering the 
correlation function 



i I d 4 xe iqx (U\Tu{x)gG ml a lb d{x) d(0) 7 „u(0)|0) = (g M g„ - g^q 2 )U C4 (q 2 ), 



which yields the following sum rule |16 



(c.i; 



^-M 2 {l-e-^ 2 (l + -^) 
18tt 3 I V M 2 J 



((a s /n)G 2 ) 
QM 2 
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27M 4 



ira s (qq)' 



(C.2) 
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In the Borel window 1 GeV 2 < M 2 < 2 GeV 2 and with so ~ 1.5 GeV 2 and the condensates 
((a s /ir)G 2 ) = (0.012 ± 0.006) GeV 4 and {^alqq) 2 = 0-56 • (-0.25 GeV) 6 , we obtain 

C 4 (/i = 1 GeV) = 0.15 ±0.10. (C.3) 

The calculation of involves dimension 6 operators for which the QCD sum rule approach 
becomes rather unreliable. Because of this, we choose to make a simple estimate by considering 
the leading contribution to a correlation function that vanishes in perturbation theory: 

CF x =i Jd 4 ye^(0\Td(y)a KX u(y)O { ^M\0) 

= - (qagGq) %{q K g XlJ - q x g KU ) + 0(l/g 4 ). (C.4) 
2 q z 

The contribution of the p meson to this correlation function is 

qu(q,gxp - qxgMX^ + V W ) + q p (q K g Xu - q x g KU )(X^ - X®) 

(C.5) 



CF X = f p fjm 3 p 
nip — q A 1 F 



from which we obtain in the local duality limit q 2 — > — oo: 

X? * -X® ^ -— ^ (qvgGq) (C.6) 



at a hadronic scale fi ~ 1 GeV. This has to be compared with the estimate for £ 4 obtained in 
the same approximation by considering a similar correlation function with the operator ( [4.18 



C < * ~2Uf^ p (C7) 
Putting in numbers, we get (4 w 0.3, which is a factor two larger than what comes from the 



more accurate (and laborious) analysis in flq| , see above. Using the definition ( p.8|) we get, 
finally 

GeV) « 7/9 (C.8) 

with, probably, a 100% error. 



C.2 Chiral-Odd 



The calculation of the matrix elements (J and (J, defined in ( 5.17Q , and the matrix elements 



((Q )) °f tne operators Q a 



(0 



1, 3, 5, defined in Q5.21 



parameters of the photon distribution function in Ref. 
paper can be adapted to the present case. 



191 



, is analogous to calculation of the 
and the sum rules obtained in this 
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To leading conformal spin accuracy, we need to estimate two parameters, Cj and (J. To 
this end, we consider the correlation functions 

CF ± = i [ d 4 ye i<iy {0\Td(y) lfl u(y)u{0)9[G{0) ± ij 5 G(0)Ud(0)\0) , (C.9) 



which vanish in perturbation theory. The leading power corrections were calculated in flS 
yielding 

CF + = i(q p g^ a - q a g^) 0(l/g 4 ), 

CF_ = i(q p g, a - q a g^) + 0(l/g 4 )| . (CIO) 

Saturation with a p meson gives, on the other hand, 

cf± = -4z~2 m %fJ(£ ± 2>(W/s - g^Pa), (c.ii) 

m p q 

so that, taking into account only the leading 1/q 2 terms, we have 

Cj + Cf = o, 

C4_C4 ""3<^- (C - 12) 

Using the same numerical input as in the last section, we obtain 

Cj(l GeV) = -Cj(l GeV) = 0.10 ± 0.05 (C.13) 

with a rather conservative large error. 

We use the same method to estimate also the ({Q^)), and consider the correlation functions 

CF^ = i J d 4 ye igy (0\Td(y) lK u(y) Q^(0)|0>. (C.14) 

As shown in [19|, the lowest order power correction to CF^ 3 ' 5 ^ vanishes, so that 

«Q(3))) = {{ Q(V)) = o (C.15) 

to that accuracy. For ((Q^)), on the other hand, the mixed condensate gives a nonzero 
contribution and we obtain 

«Q<'>>>(1 GeV) ^<^@^. -0.30 (C.16) 

J m pJpJp 

This value is likely to be overestimated since the mass scale in the correlation function is much 
larger than the p meson mass, see the discussion in |f[9 |. We thus prefer to give 

((Q (1) »(l GeV) = -0.15 ±0.15 (C.17) 

clS db conservative estimate. 
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